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We use neutron resonance spin echo and Larmor diffraction to study the effect of uniaxial pres¬ 
sure on the tetragonal-to-orthorhombic structural (T s ) and antiferromagnetic (AF) phase transitions 
in iron pnictides BaFe 2 -a;Ni x As 2 (x — 0,0.03,0.12), SrFei. 9 7 Nio.o 3 As 2 , and BaFe 2 (Aso. 7 Po. 3 ) 2 - In 
antiferromagnetically ordered BaFe 2 - cc Ni ;r As 2 and SrFei. 97 Nio.o 3 As 2 with Tn and T s (Tn < T s ), a 
uniaxial pressure necessary to detwin the sample also increases Tn , smears out the structural transi¬ 
tion, and induces an orthorhombic lattice distortion at all temperatures. By comparing temperature 
and doping dependence of the pressure induced lattice parameter changes with the elastoresistance 
and nematic susceptibility obtained from transport and ultrasonic measurements, we conclude that 
the in-plane resistivity anisotropy found in the paramagnetic state of electron underdoped iron 
pnictides depends sensitively on the nature of the magnetic phase transition and a strong coupling 
between the uniaxial pressure induced lattice distortion and electronic nematic susceptibility. 

PACS numbers: 74.70.Xa, 75.30.Gw, 78.70.Nx 


I. INTRODUCTION 

The parent compounds of iron pnictide superconduc¬ 
tors such as BaFe2 As 2 and SrFe2 As2 exhibit a tetragonal- 
to-orthorhombic structural transition at T s followed by 
development of collinear antiferromagnetic (AF) order 
along the a-axis of the orthorhombic lattice below T/v 
[left inset in Fig. 1 (a) and T s » TvpHl. Upon electron¬ 
doping via partially substituting Fe by Co or Ni to form 
BaFe2- rE T[ c As2 (T = Co, Ni), the nearly coupled struc¬ 
tural and magnetic phase transitions in BaFe2As2 be¬ 
come two separate second order phase transitions at 
T s and T/v (T s > Tn) that decrease in temperature 
with increasing x 7 10 . On the other hand, the coupled 
first order structural and magnetic phase transitions in 
SrFe2As2 4 [ while decreasing in temperature with increas¬ 
ing x in SrFe2- cc T c As2, remain coupled first order tran¬ 
sitions leading up to superconductivity 11 . 

Because the structural and magnetic phase transi¬ 
tions in BaFe2- cc T c As2 and SrFe2- cc T c As2 occur below 
room temperature, iron pnictides in the orthorhombic 
AF ground state will form twin domains with AF Bragg 
peaks appearing at the in-plane (± 1 , 0 ) and ( 0 , ± 1 ) po¬ 
sitions in reciprocal space [right inset in Fig. 1 (a)] 6 . 
To probe the intrinsic electronic properties of these ma¬ 
terials, one can apply uniaxial pressure along one axis 
of the orthorhombic lattice to obtain single domain 
samples 12 15 . Indeed, transport measurements on uni¬ 
axial pressure detwinned electron-doped BaFe2- x T x As2 
(T = Co, Ni) reveal in-plane resistivity anisotropy in the 
AF state that persists to temperatures above the zero- 


pressure T/v and T s 12 ®^. On the other hand, similar 
transport measurements on uniaxial pressured detwinned 
SrFe2-xTcAs2 (T = Co, Ni) indicate vanishingly small 
resistivity anisotropy at temperatures above the zero 
pressure coupled T/v and Figure 1 (b) compares 

temperature dependence of the resistivity anisotropy [de¬ 
fined as A p = (p b - Pa)/(Pb + Pa), where p a and p b are 
resistivity along the a and b axis of the orthorhombic lat¬ 
tice, respectively] obtained under 20 MPa uniaxial pres¬ 
sure for BaFe 2 As 2 , BaFei. 9 7 Nio.o 3 As 2 , SrFei. 9 7 Nio.o 3 As 2 , 
and SrFe2As2- Consistent with earlier works lJHHI we find 
that resistivity anisotropy is much larger in BaFe2As2 
and BaFei.g 7 Nio.o3As2 at temperatures above T/v- 

Although resistivity anisotropy in the paramagnetic 
state of the iron pnictides under applied uniaxial pres¬ 
sure suggests the presence of an electronic nematic phase 
that breaks the in-plane fourfold rotational symmetry 
(C4) of the underlying tetragonal latticd^HIS much is 
unclear about the microscopic origin of the in-plane re¬ 
sistivity anisotropy and electronic nematic phas^SHSHl^ 
Since neutron scattering experiments reveal that uni¬ 
axial pressure necessary to detwin the sample also in¬ 
creases T/v of the system, the observed in-plane resistivity 
anisotropy above the zero pressure T/v and T s may arise 
from the increased Tn and intrinsic anisotropic nature 
of the collinear AF phas^ 39 ®. Furthermore, while it is 
generally assumed that the uniaxial pressure for sample 
detwinning has negligible effect on the lattice parame¬ 
ters of the iron pnictides 12 ^, the precise effect of uni¬ 
axial pressure on structural distortion of these materi¬ 
als is unknown. From neutron extinction effect measure- 
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merits, a uniaxial pressure is suggested to push structural 
fluctuations related to the orthorhombic distortion to a 
temperature well above the zero-pressure value of T<pl, 
similar to the effect on the resistivity anisotropy 12 15 . 
To understand the microscopic origin of the in-plane 
resistivity anisotropy in the paramagnetic state^HED it 
is important to establish the effect of a uniaxial pres¬ 
sure on the magnetic and structural phase transitions 
of BaFe2-^T E As2 and SrFe2-^T E As2, and determine if 
the electronic anisotropy in the pa ramag netic tetragonal 
phase of iron pnictides is intrinsic!- 42 * 43 ^ , or entirely due 
to the symmetry breaking uniaxial pressure applied to 
the material d 44 * 45 *. It is also important to deduce what 
role the natur e of th e AF transition plays in the nematic 
susceptibiht}^ ^ 26 * 38 * and how the latter depends on the 
uniaxial pressure. 


In this paper, we use neutron resonance spin echo 
(NRSE^pSEl and Larmor diffractiorP^ to study the 
effect of uniaxial pressure on the structural and 
magnetic phase transitions in electron doped iron 
pnictides BaFe2- ;E NL c As2 with x = 0 , 0 . 03 , 0 .l^P^ 

and SrFei.97Nio.Q3AsJ^j, and in the isovalently 
doped BaFe2(Aso.7Po.3)2r^- While the underdoped 
BaFei.97Ni 0 .o3As2 (TJv = 109 K and T s = 114 K) exhibits 
a second-order AF transition below T s , SrFei.97Nio.o3As2 
has coupled first-order structural and magnetic phase 
transitions at TV = T s « 175 KP^. The electron over¬ 
doped BaFe1.88Nio.12A.s2 (T c = 18.6 K) and isovalently 
doped BaFe2(Aso.7Po.3)2 (T c = 30 K) have a paramag¬ 
netic tetragonal structure at all temperatures without 
static AF order. Figure 1 (c) summarizes the key experi¬ 
mental result of the present work, where the temperature 
dependences of the uniaxial pressure induced orthorhom¬ 
bic lattice distortion S(P = 20 MPa) — 5 (P = 0 MPa) 
are determined using neutron Larmor diffraction for 
BaFe 2 As 2 , BaFei.97^0.03 As 2 , SrFei.9 7 Ni 0 .o3As2, and 
BaFei.88Nio.i2As2 [we defined the lattice distortion 
S = (a — b) / (a + b) with a and b being the orthorhom¬ 
bic lattice parameters]. Remarkably, the magnitude 
of our determined structural nematic susceptibility 
dS/dP oc S(P) — 5 ( 0 ) in Figure 1 (c) is comparable 
in all three materials that have a structural phase 
transition, unlike the very different values of the resis¬ 
tivity anisotropy displayed in Figure 1 (b). Comparing 
these results with those of the elastoresistan ce an d 
nematic susceptibility obtained from transport 20 23 ^ 
and from elas tic shea r modulus/ultrasound spectroscopy 
measurement ^M 15 i| 52 j, we conclude that the resistivity 
anisotropy in the paramagnetic phase of the iron pnic¬ 
tides depends sensitively on whether the underlying 
magnetic phase transition is first or second order. We 
also find a strong coupling between the uniaxial pressure 
induced lattice distortion and the electronic nematic 
susceptibility, and have to be cautious in directly relating 
resistivity anistropy to the nematic order parameter in 
the iron pnictides. 



Temperature (K) 

FIG. 1: (Color online) (a) The schematic electronic phase dia¬ 
gram of BaFe 2 -a;NLAs 2 with arrows marking x = 0, 0.03 and 
0.12 samples described in the present study. The AF, PM, 
Ort, Tet, IC, SC are antiferromagnetic, paramagnetic, or¬ 
thorhombic, tetragonaLincommensurate, and superconduct¬ 
ing states, respectively 10 . The left inset shows the direction 
of the applied uniaxial pressure (marked by the vertical ar¬ 
rows) and the spin arrangements of Fe in the AF ordered 
iron pnictides, where a and b are the orthorhombic axes. The 
right inset shows the corresponding reciprocal lattice. All 
the marked positions have AF or nuclear Bragg peaks for a 
twinned sample, while the positions marked by open sym¬ 
bols have vanishing scattering intensity for a detwinned sam¬ 
ple. (b) Temperature dependence of the resistivity anisotropy 
for BaFe 2 -ccNi x As 2 and SrFe 2 -xNi x As 2 {x — 0,0.03) under 
P « 20 MPa. (c) Summary of temperature dependence of 
the uniaxial pressure induced lattice distortion at P = 20 
MPa [S(P = 20 MPa) — S(P = 0 MPa)] for BaFe 2 - ;r Ni ;r As 2 
(x — 0,0.03,0.12) and SrFei. 97 Nio.o 3 As 2 . The actual data for 
x — 0.03, 0.12 are normalized to 20 PMa assuming a linear re¬ 
lationship between uniaxial pressure and 5. Uniaxial pressure 
induced lattice distortion vanishes rapidly below Tn marked 
by the vertical dashed lines in (b) and (c). 


II. RESULTS 
A. Experimental Results 

Our experiments were carried out using conventional 
thermal triple-axis spectrometer PUMA and three axes 
spin echo spectrometer (TRISP) at the Forschungsneu- 
tronenquelle Heinz Maier-Leibnitz (MLZ), Garching, 
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FIG. 2: (Color online) (a) Magnetic order parameters at 
Q = (1,0,1) for the zero (P = 0) and uniaxial pressured 
(P ~ 15MPa) BaFei. 97 Nio.o 3 As 2 - Tjv is 109K for an unpres¬ 
sured sample (blue diamonds). Upon applying uniaxial pres¬ 
sure of P ~ 15 MPa, the Tn is enhanced to 118K and the sam¬ 
ple becomes 100% detwinned as seen by PUMA and TRISP 
measurements, (b) The energy line width (Half-Width-at- 
Half-Maximum, T) of the magnetic Bragg peak Q = (1,0,1) 
measured by NRSE using TRISP for BaFei. 97 Nio.o 3 As 2 . The 
blue and red dashed lines indicate Tn in P = 0 and 15 MPa 
unaixial pressure, respectively. The slight larger errors of T 
near Tn is due to low statistics data. 


Germany. The principles of NRSE and Larmor diffrac¬ 
tion are described elsewhere 50 . Single crystals of 
BaFe 2 - cc Ni a ,As 2 , SrFe 2 _ x Ni x As 2 , and BaFe 2 (As 0 . 7 Po 3)2 
were grown by self-flux method as described beford^^. 
We define the momentum transfer Q in the three- 
dimensional reciprocal space in A -1 as Q = ida*+Pb* + 
Pc*, where id, iU, and L are Miller indices and a* = 
a27r/a, b* = b27r/6, c* = c2tt/c with a « b « 5.6 A, and 
c = 12.96 A for BaFe 2 _ a ,Ni ;E As 2 . In this notation, the AF 
Bragg peaks should occur at (±1, 0, L ) (P = 1, 3, 5, • • •) 
positions in reciprocal space of a completely detwinned 
sample [right inset in Fig. 1(a)]. For neutron scatter¬ 
ing experiments, single crystals are aligned in either the 
[id, K,H + iFp or [id, K, 0] zone. 

We first discuss the effect of uniaxial pressure on the 
collinear AF order in BaFe 2 _ ;E T ;E As 2 . In previous neu¬ 
tron scattering work on BaFe 2 _ a ,Co ;E As 2 , the Neel tem¬ 
perature (T/v) was found to be pushed to higher tem¬ 
perature under uni axia l strain field, forming a broader 
magnetic transitior i 39 * 4Q [ Moreover, it seems that t he in - 
crease in T/v depends on the annealing condition 41 54 . 
Although the T/v enhancement was attributed to uniax¬ 
ial strain aligned fluctuating magnetic domains, the effect 
of uniaxial pressure on the ordered moment remains elu¬ 
sive and the nature of the T/v enhancement is still under 
debat e™^ . 

By aligning single crystals in the [1, 0,1] x [0,1,1] scat¬ 
tering planep we were able to determine T/v, detwinning 
ratio, as well as the ordered moment of the system un¬ 


der zero and finite uniaxial pressures. Figure 2(a) shows 
temperature dependence of the (1, 0,1) and (0,1,1) mag¬ 
netic scattering intensity for BaFei.g 7 Nio.o 3 As 2 obtained 
using PUMA [left axis in Fig. 2(a)] and TRISP (right 
axis). The two sets of data are in excellent quantitative 
agreement with each other. Under the applied uniaxial 
pressure of P « 15 MPa, the Neel temperature increases 
from T/v ~ 109 K (at P = 0) to Tn ~ 118 K. The mag¬ 
netic scattering intensity [Fig. 2(a)] in the (1,0,1) peak 
becomes approximately twice as large as in the twinned 
sample, whereas the (0,1,1) peak vanishes, suggesting 
that the sample is completely detwinned and the applied 
uniaxial pressure does not significantly affect the ordered 
moment. 

To test whether the T/v increase is an intrinsic fea¬ 
ture of the system, we note that the magnetic order pa¬ 
rameter under uniaxial pressure has a round tail around 
T/v I 39 I 4 Q I suggesting that the T/v enhancement could arise 
from enhanced slow spin dynamics (critical scattering) 
under inhomogeneous uniaxial strain field and cannot be 
resolved by conventional triple-axis neutron diffraction 
due to its coarse energy resolution (A E ~ 0.3 — 1 meV). 
To clarify the nature of the increase in T/v, we have mea¬ 
sured the energy line-width (T > 0, [see Fig. 2(b)]) of 
the quasielastic scattering for magnetic reflection (1,0,1) 
using hig h energy resolution (A E ~ 1 /ieV) NRSE at 
TRISP^. As seen in Figure 2(b), the T at all measured 
temperatures are resolution limited, indicating that the 
increase in magnetic scattering intensity below Tn ~ 118 
K is elastic (T < 1 /ieV), and an intrinsic nature of the 
system. 

To determine the effect of uniaxial pressure on the 
tetragonal-to-orthorhombic phase transition in iron pnic- 
tides, we carried out neutron Larmor diffraction experi¬ 
ments capable of measuring minor change of lattice spac¬ 
ing d = 27r/|Q(iJ, id, L)\ and its spread Ad with a reso¬ 
lution better than 10 -5 in A d/d [inset in Fig. 3(a) ] 48 | 5Q | . 
We focus on (4, 0, 0) and (0,4, 0) nuclear Bragg reflections 
corresponding to a d-spacing d ~ a/4, which we measured 
in BaFe 2 _ a ,NL c As 2 (x = 0,0.03,0.12), SrFei. 9 7 Ni 0 .o 3 As 2 , 
and BaFe 2 (Aso. 7 Po. 3)2 both on freshly prepared samples 
(uniaxial pressure P = 0) and under uniaxial pressure 
(P « 10,20 MPa). Figure 3 shows the temperature 
and pressure dependence of the d spread for these sam¬ 
ples. The d spread are characterized by the FWHM 
(Full-Width-Half-Maximum) of the lattice spacing dis¬ 
tribution /(Ad/d), which is assumed to be Gaussian 
distribution^ 0 !. The diamonds in Figure 3(a) show tem¬ 
perature dependence of the FWHM for BaFe 2 As 2 at zero 
pressure. Similar to BaFei.g 7 Nio.o 3 As 2 p^, temperature 
dependence of FWHM follows a Curie-Wiess form and 
peaks around the zero-pressure value of T/v ~ T s . Upon 
application of a uniaxial pressure P ~ 20 MPa, the mag¬ 
nitude of FWHM increases at all temperatures and now 
peaks at an enhanced T/v = 144 K [Fig. 3(a)]. 

Figure 3(b) shows similar data for SrFei.g 7 Nio.o 3 As 2 , 
where there are coupled strong first order structural and 
AF phase transitions at T/v = T s = 175 kPA Compared 
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FIG. 3: (Color online) Temperature dependence of FWHM 
of A d/d in sev eral iron pnictides under different uni¬ 
axial pressures obtained from neutron Larmor diffraction 
experiments 50 ! (a) Temperature dependence of A d/d in 
FWHM for the (4, 0, 0) Bragg reflection of BaFe 2 As 2 at P — 0 
and 20 MPa. The solid line above Tn is a fit using Curie-Wiess 
formalism [FWHM(T) = A/(T — Ti) + P, where A, B, and T\ 
are fitting parameters], (b) Similar data for SrFei.g 7 Nio.o 3 As 2 . 
The vertical blue and red dashed lines in (a) and (b) mark 
the Tn of the sample at zero and finite pressure, respec¬ 
tively. (c) Similar data for BaFei. 9 7 Nio.o 3 As 2 , where the ver¬ 
tical green and blue dashed lines mark Tn and T s , respec¬ 
tively, at zero pressure. The open green and yellow squares 
mark measurements of FWHM under zero pressure (fresh) 
at the (4,0,0) and (0,4,0) Bragg peaks, respectively. The 
pink solid circles are identical measurements under P « 10 
MPa uniaxial pressure on (4, 0, 0). The solid green diamonds 
are data after uniaxial pressure is released. The vertical red 
dashed line indicate the peak position of the FWHM under 
P — 10 MPa. (d) Temperature dependence of FWHM in 
A d/d for BaFe 2 (Aso. 7 Po. 3)2 at P = 0 MPa (solid green cir¬ 
cles), BaFei. 88 Nio.i 2 As 2 at P = 0 (solid green diamonds) and 
14 MPa (solid red circles). 


with BaFe 2 As 2 , where the AF phase transition is weakly 
first order and structural transition is second order^, the 
AF and structural transition induced changes in FWHM 
are much smaller and confined to temperatures close to 
T n « T s in SrFei. 9 7 Nio.o 3 As 2 [Fig. 3(b)]. Under a uni¬ 
axial pressure P « 20 MPa, however, both the FWHM 
and T/v increase dramatically with solid lines showing 
Curie-Wiess fits to the data. For BaFei. 97 Nio.osAs 2 , ap¬ 
plication of a P ^ 10 MPa uniaxial pressure transforms 
temperature dependence of the FWHM, which forms a 
broad peak above the zero-pressure value of T s . Upon 
releasing the uniaxial pressure [P released, filled green 
diamonds in Fig. 3(c)], the system goes back to the orig¬ 
inal unpressured fresh state. 

Figure 3(d) compares temperature dependence of the 
FWHM for electron overdoped BaFei. 8 sNio.i 2 As 2 and 
BaFe 2 (Aso. 7 Po. 3 ) 2 , where both materials are in the para¬ 
magnetic tetragonal state without static AF order. The 
weak temperature dependence of FHWM in these ma¬ 
terials suggests that the large temperature dependence 
of FWHM in AF ordered BaFe 2 -^NUAs 2 (x = 0, 0.03) 
and SrFei. 97 Nio.o 3 As 2 is due to a strong magnetoelastic 
coupling. Although application of a P « 14 MPa uniax¬ 
ial pressure on BaFei. 8 sNio.i 2 As 2 increases the absolute 
value of FWHM, it is still weakly temperature dependent 
[Fig. 3(d)]. 

To further demonstrate the impact of uniaxial pres¬ 
sure on the tetragonal-to-orthorhombic structural tran¬ 
sition in BaFe 2 - :E NL c As 2 (x = 0,0.03,0.12) and 

SrFei. 9 7 Ni 0 .o 3 As 2 , we compare in Figure 4 tempera¬ 
ture dependence of the lattice parameters along the 
orthorhombic a and b axis directions under zero and 
finite uniaxial pressure. We first discuss results for 
BaFe 2 -a;NUAs 2 with x = 0 [Fig. 4(a), 4(b)] and 0.03 
[Fig. 4(c) and 4(d)]. At P = 0, the lattice parame¬ 
ters have a = b at temperatures above T s (tetragonal 
phase) and decrease linearly with decreasing tempera¬ 
ture [open diamonds and hexagons in Fig. 4(a) and 
4(c)]. Upon application of a uniaxial pressure, the sys¬ 
tem becomes orthorhombic at all temperatures and the 
orthorhombic structural transition becomes a crossover 
[filled diamonds and hexagons in Fig. 4(a) and 4(c)]. 
Figures 4(b) and 4(d) show temperature dependence of 
the lattice orthorhombicity 5 = (a — b)/(a + b) at differ¬ 
ent uniaxial pressures for x = 0, and 0.03, respectively. 
For unpressured fresh samples (P = 0), and after the 
pressure has been released, the tetragonal structure be¬ 
comes orthorhombic below T s and the AF order below 
T/v further enhances the lattice orthorhombicit}®. Upon 
applying the uniaxial pressure P « 10, 15, and 20 MPa, 
the temperature dependence of the lattice orthorhombic¬ 
ity becomes remarkably similar to that of the B 2 g elas- 
toresistance and nematic susce ptibility of BaFe 2 - x T a; As 2 
obtained from transport^ 3 ^ and elastic shear m od- 
ulus/ultrasound spectroscopy measurement^^^, re¬ 
spectively. 
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B. Theoretical Ginzburg—Landau analysis 


To understand the temperature dependence of the 
pressure-induced lattice orthorhombicity described in 
Figs. 4(b) and 4(d), we consider the Ginzburg-Landau 
free energy formalism used in previous works 20 38 : 

F[p, 8 ) * F 0 +1 (T - To )- A 5<p - PS, 

( 1 ) 

where the electronic nematic order parameter p is cou¬ 
pled linearly to the orthorhombic lattice distortion 5. It 
then follows that (see appendbP^) 


{=(A^)+P)/C 66)0j (2) 


where Cq e,o is the bare elastic constant that has no strong 
temperature dependence and P is the conjugate uniaxial 
pressure (stress )^ 8 l 41 l 5Q t t52] j n ^h.e absence of the elasto- 
nematic coupling (A = 0), the nematic susceptibility 
Xcp = l/[a(T — To)] is characterized by the Curie-Weiss 
temperature Tq. Upon considering the coupling between 
the nematic order parameter p and the structural lattice 
distortion S (or equivalently, the elastic s hear s train £ 6 ), 
the elastic susceptibility takes on the form 20 38 : 


d 5 _ 1 T-T 0 

d P ~ C^T-T^w 


( 3 ) 


with the renormalized nematic transition temperature 
Tf w = T 0 + A 2 /(aC 66) o) that is increased compared 
to the bare Curie-Weiss temperature Tq. The pressure- 
induced lattice distortions in Figs. 4(b) and 4(d) can 
be well described by the Curie-Weiss functional form 50 . 
Therefore, uniaxial pressure induced orthorhombic lat¬ 
tice distortion and its temperature dependence in un¬ 
doped and underdoped BaFe 2 - 2C Ni ;E As 2 are directly as¬ 
sociated with the nematic susceptibilit y 24 1 38 1 . Since the 
external uniaxial pressure explicitly breaks the tetrago¬ 
nal lattice symmetry, it turns the nematic transition at 
T s = T^ w into a crossover, as is clearly seen in Fig. 4. 

If the in-plane resistivity anisotropy in electron un¬ 
derdoped iron pnictides indeed arises from the cou¬ 
pling of the uniaxial-pressure induced lattice distortion 
5 with the nematic susceptibility, it would be interest¬ 
ing to determine the effect of similar uniaxial pressure 
on the electron overdoped sample, where the resistiv¬ 
ity anisotropy is known to be much weakei^. Fig¬ 
ures 4(e) and 4(f) summarize the outcome of the neu¬ 
tron Larmor diffraction experiments on uniaxial pres¬ 
sured BaFe 1 . 88 Nio. 12 A.s 2 , which is tetragonal (a = b ) and 
non-magnetic at all temperatures in zero pressurd^. Fig¬ 
ure 4(e) shows temperature dependence of the lattice 
parameter changes along the a-axis (A a/a) and 6-axis 
(A b/b) under a uniaxial pressure of P ~ 14 MPa. For 
comparison, the thermal contraction of aluminum is also 
shown®!. Figure 4(f) shows the temperature dependence 
of the orthorhombic lattice distortion S, which reveals 
a clear anomaly at T c co nsistent with ultrasonic spec¬ 
troscopy measurements 51 52 . While the applied uniaxial 


pressure induces orthorhombic lattice distortion at 230 
K, the magnitude of the lattice distortion, S « 1.1 x 10 -4 , 
is about 5 times smaller than that of BaFe 2 As 2 and 
BaFei. 97 Nio.o 3 As 2 at 230 K. On cooling to 20 K, S in 
BaFei. 88 Nio.i 2 As 2 increases to ^ 2 x 10 -4 , while S in 
BaFe 2 As 2 and BaFei. 97 Ni 0 . 03 As 2 becomes ^ 2.5 x 10 -3 
near T s [Fig. 4(b) and 4(d)], an order of magnitude larger 
than that of the electron overdoped compound. 

To understand how a uniaxial pressure affects the first 
order nature of the structural and magnetic phase tran¬ 
sitions in SrFei .97Nio. 03 AS 2 , we compare in Fig. 4(g) and 
4(h) temperature dependence of the lattice parameters 
and orthorhombicity under the zero and finite uniax¬ 
ial pressure. At zero pressure, the first order nature of 
the structural transition is clearly seen in hysteresis of 
temperature dependence of the lattice parameters and 
distortion [Fig. 4(g) and 4(h)]. Upon application of 
P ~ 20 MPa uniaxial pressure, the lattice orthorhom¬ 
bicity no longer displays the first order transition at 
T s , but instead becomes a crossover, similar to that ob¬ 
served in the undoped and underdoped BaFe 2 -ccNi x As 2 
[see Figs. 4(b) and 4(d)]. 

Assuming that the application of the modest uniax¬ 
ial pressure P « 20 MPa can be considered in the 
linear-response regime®!, we can estimate the elastic 
susceptibility from the finite difference dS/dP oc A(S) 
= 6 (P = 20 MPa) — S(P = 0 ) and compare it among 
the different compounds in the iron pnictide fam¬ 
ily. Figure 1(c) compares temperature dependence of 
5(P = 20 MPa) — S(P = 0) for BaFe 2 - ;E NL c As 2 (x = 
0,0.03,0.12) and SrFei. 97 Nio.o 3 As 2 normalized for P = 
20 MPa. For AF ordered BaFe 2 -a;Ni x As 2 (x = 0, 0.03) 
and SrFei. 97 N 10 . 03 AS 2 , the magnitudes of the pressure- 
induced lattice orthorhombicity are similar in the para¬ 
magnetic phase and vanish rapidly upon entering into the 
AF ordered state. Furthermore, the S(P = 20 MPa) — 
S(P = 0 ) decreases for the iron pnictides with reduced 
T/v, and are much smaller for BaFei. 8 sNio.i 2 As 2 . 


III. DISCUSSION 

It is well known that the effect of increasing electron¬ 
doping in BaFe 2 - cc T E As 2 is to suppress the static AF 
order and to eliminate the low-temperature lattice 
orthorhombicity 7 A3. At zero pressure, BaFe 2 As 2 first 
exhibits a second-order structural transition from the 
high-temperature paramagnetic tetragonal phase to a 
paramagnetic orthorhombic phase at T s , followed by 
a discontinuous further orthorhombic structural distor¬ 
tion and weakly first order AF phase transition at T/v 
(T/v < T s ) due to magnetoelastic coupling 5 . Upon Ni- 
doping in BaFe 2 - ;E NL c As 2 , the structural and magnetic 
phase transitions are gradually separated and suppressed 
[Fig. 1(a)], and become second order in natur^HSI. 
Upon application of a uniaxial pressure, the C 4 rota¬ 
tional symmetry of the tetragonal lattice is broken. Since 
the tetragonal-to-orthorhombic symmetry of the under- 
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FIG. 4: (Color online) Temperature dependence of the a and b lattice parameters and orthorhombicity 5 under different 
uniaxial pressure conditions (P = 0 fresh, ~10, ~15, ~20 and 0 released MPa) for BaFe 2 -^Ni cc As 2 (x = 0,0.03,0.12) and 
SrFei. 97 Nio.o 3 As 2 . (a) Temperature dependence of the a and b under P — 0 and 20 MPa uniaxial pressure for BaFe 2 As 2 . 
(b) Temperature dependence of 6 under different uniaxial pressure (P = 0 fresh, ~20, and 0 released MPa). The vertical 
blue dashed line marks the TJv /T s . (c) Temperature dependence of the lattice parameters a and b at P = 0 and 15 MPa for 
BaFei. 97 Ni 0 . 03 As 2 . (d) Temperature dependence of S under different uniaxial pressure (P = 0 fresh, ~10, ~15 and 0 released 
MPa). The open red diamonds and green hexagons are obtained by lattice thermal expansion measurements under uniaxial 
pressure. The pink squares are measurements of an unpressured fresh sample, while blue circles are obtained after releasing 
P w 10 MPa uniaxial pressure. The blue circles and two pink squares below T s are from zero pressure Larmor diffraction 
measurements. The pink squares above T s are obtained from thermal expansion measurements. The vertical green and blue 
dashed lines in (c) and (d) mark the T/v and T s of the sample at zero pressure, respectively, (e) Temperature dependence of 
the a and b for BaFei. 8 sNio.i 2 As 2 . The lattice thermal expansion of aluminum is plotted as a reference 56 . The vertical arrow 
marks T c = 18.6 K and the solid lines are guides to the eye. (f) Temperature dependence of the orthorhombic lattice distortions 
calculated from (c). (g) Temperature dependence of the a and b lattice parameters for SrFei. 97 Nio.o 3 As 2 under P = 0 fresh, 20 
and 0 (released) MPa. (h) Temperature dependence of 5 for the same pressure condition. The solid curves in (b), (d), (f), and 
(h) are fits using a Curie-Weiss functional forrrP^. 


lying lattice can only be broken once, T s will become a 
crossover regardless the magnitude of the applied pres¬ 
sure, as our findings in Figs. 3 and 4 corroborate. The 
same conclusion holds for SrFei. 97 Nio.o 3 As 2 where the 
first-order structural transition becomes a crossover [see 
Figs. 4(g) and 4(h)]. Therefore, both BaFe 2 -zNT c As 2 


and SrFei. 97 Nio.o 3 As 2 under uniaxial pressure can only 
exhibit AF phase transition. We note that our mea¬ 
surements and theoretical Landau-Ginzburg analysis do 
not rely on the microscopic nature of the nematic order 
parameter p. In particular, they apply equa lly well to 
the so-called Ising spin nematic scenario 25 29 or the or- 
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bital order interpretation of nematicit}®®. In fact, the 
ferro-or bital order ip^ = (n xz — n yz ) is always linearly 
couplecPSIM] to the Ising spin nematic order parameter 
c Pspin = (Si • S i+z — Si • Si+y), so the orbital order is 
generically present whenever (p sp i n 7 ^ 0 , although there 
are theoretical indications that the converse is not always 
true. In other words, the orbital order can exist in the 
absence of static AF orde]®} as is known to be the case in 
FeSe 57 59 . In either case, the application of external uni¬ 
axial stress renders the nematic transition a crossover, so 
that the lattice distortion S and consequently ip are both 
finite above the zero-stress value of T s . In this light, the 
electronic anisotropy seen in the magnetic torque 21 and 
scanning tunneling microscop}^ measurements above T s 
without explicit external uniaxial pressure is likely due to 
intrinsic local strain in these materials which breaks the 
C 4 rotational symmetry of the paramagnetic tetragonal 
phase. Indeed, local strain-induced effect has recently 
been observed in free standing BaFe 2 As 2 above T/v and 


The key finding of the present work is that un¬ 
doped BaFe 2 As 2 , as well as BaFei. 97 Nio.o 3 As 2 and 
SrFei. 97 Nio.o 3 As 2 all exhibit similar magnitudes of the 
pressure-induced lattice orthorhombicity [Fig. 1(c) and 
Fig. 4] and FWHM of A d/d near T/v (Fig. 3). This in¬ 
dicates that these samples experience similar strain field 
under nominally similar applied uniaxial pressure, thus 
suggesting that the doped Ni impurities do not play an 
important role in determining the strain field inside the 
sample. Theoretically, the electronic anisotropy of the 
iron pnictides is expected to couple linearly to the lat¬ 
tice orthorhombicity as ca ptured by the effec¬ 

tive Landau free energy in Eq. §• The Curie-Weiss 
like temperature dependence of the uniaxial pressure in¬ 
duced lattice distortion [Fig. 1(c)] is consistent with 
the temperature dependence of the nematic suscepti¬ 
bility dS/dP in Eq. (3) and agrees with the results of 
Young’s modulus measurementThis gives us con¬ 
fidence that in the effective Landau description^!, the 
uniaxial pressure-induced lattice distortion S has a com¬ 
ponent proportional to the electronic nematic order pa¬ 
rameter ip via Eq. ( 2 ), where one expects S oc ip in zero 
pressure (P = 0 ). Since S has similar magnitude in 
BaFe 2 As 2 , BaFei. 97 Nio.o 3 As 2 , and SrFei. 97 Nio.o 3 As 2 (see 
Fig. 4), one would also expect comparable values of ip 
in all three compounds. So if one uses the resistivity 
anisotropy A p = (p^ — p a )/{pb + Pa) as a proxy for the 
nematic order parameter, as has been widely used in the 
literature 12 2 Q I 23 I 24 ( how does one then explain the resis¬ 
tivity anisotropy differences in BaFc 2 - T x As 2 14 and a 
much smaller resistivity anisotropy above T/v [Fig. 1 (b)] 
in SrFe 2 -^TcAs 2 family of materiald^? The bare value 
of the elastic shear modulus that enters Eq. ( |20| ) 

has no strong temperature dependenc ^ 8 and from the 
Curie-Weiss fits of the nematic susceptibility to Eq. (24), 
we find it to be roughly the same in all three com¬ 
pounds, C 66,0 ~ 50 GPaP^. The only remaining un¬ 
known variable is the elasto-nematic coupling constant A, 


which could be material-dependent but not temperature- 
dependent 20 38 . It is thus very challenging to explain 
the qualitatively different temperature dependence of the 
resistivity anisotropy in BaFei. 97 Nio.o 3 As 2 [monotonic, 
blue diamonds in Fig. 1(b)] from that in BaFe 2 As 2 and 
in SrFei. 97 Nio.o 3 As 2 [both non-monotonic, with a maxi¬ 
mum at or just below T s \. One possible explanation for 
the non-monotonic temperature dependence of the re¬ 
sistivity anisotropy, recently proposed in the context of 
FeSe®} is to assume a temperature dependent coefficient 
of proportionality between A p and ip: 

A Pin = T(TMT), (4) 

such that Y(T) tends to zero as T 0, whereas p(T) 
is expected to increase monotonically below T s as the 
temperature is lowered (consider for instance the mean- 
field result p(T) oc y/T s — T for the second order phase 
transition). 

Even with the introduction of Y(T) in Eq. , which 
has a meaning of the temperature-dependent scatter¬ 
ing function, it is extremely difficult to explain the 
much lower value of A p in SrFei. 97 Nio.o 3 As 2 compared 
to BaFe 2 As 2 and BaFei. 97 Nio.o 3 As 2 . In fact, from the 
Curie-Weiss fits of the susceptibility data, we estimate 
the elasto-nematic coupling constant A to be a fac¬ 
tor of ^ 5 smaller in SrFei. 97 Nio.o 3 As 2 compared to 
BaFei. 97 Nio.o 3 As 2 p^. Given the comparable magnitudes 
of 5 between the two compounds [see Fig. 4(d) and 
4(h)], one would then expect the nematic order parame¬ 
ter ip to be a factor of ^ 5 greater in SrFei. 97 Nio.o 3 As 2 , 
to ensure that the left-hand side of Eq. (2) remains of 
the same magnitude. And yet the resistivity anisotropy 
A p oc ip paints a diametrically opposite picture, being 
much smaller in SrFei. 97 Nio.o 3 As 2 - 

We propose that a likely resolution of this dilemma 
lies in the nature of the magnetic phase transition which 
we have so far neglected in our analysis. Indeed, it 
is well established that structural and magnetic phase 
transitions in SvFe 2 - x T x As 2 are coupled first order tran¬ 
sitions that decrease with increasing x before vanish¬ 
ing near optimal superconductivity^} while electron- 
doped BaFe 2 As 2 has second order magnetic and struc¬ 
tural phase transitions 7 10 . Although application of 
a uniaxial pressure renders the structural transition a 
crossover, the first order nature of the magnetic transi¬ 
tion means a vanishing critical regime with suppressed 
low-energy spin fluctuations at temperatures near T/v, 
compared with those of BaFe 2 - a ,T E As 2 where the AF 
phase transition is second order. One expects the scat¬ 
tering of electrons on the magnetic fluctuations, and 
hence the resistivity, to therefore be smaller in the vicin¬ 
ity of the first-order magnetic transition, as is the case 
in SrFe 2 - cc T c As 2 . We thus conclude that the vanish¬ 
ing resistivity anisotropy above T/v in the uniaxial pres¬ 
sure detwinned SrFe 2 - a ,T E As 2 (compared with those of 
BaFe 2 - a ,T E As 2 ) is likely rooted in the first order na¬ 
ture of the AF phase transition. This is also consistent 
with the increased paramagnetic resistivity anisotropy on 










moving from BaFe 2 As 2 to BaFe 1 . 96 Coo .04 a S 2 R where 
the magnetic transition changes from weakly first or¬ 
der to second ordei®^. Similarly, the lack of large 
resistivity anisotropy in the paramagnetic state of uni¬ 
axial pressured Bai_ a; K, E Fe 2 AsJ^l, Bai_ a ,Na ;E Fe 2 As^, 
and Cai_ ;r La x Fe 2 As 2 f^is likely due to the first order na¬ 
ture of the paramagnetic to AF phase transition in these 
materials. The phenomenological Landau theory can be 
extended to include the coupling betw een nematicity cp 
and the magnetic order parameter 38 * 50 1 , and our theoreti¬ 
cal analysis showd^ that the resulting uniaxial pressure- 
induced lattice distortion S(P) — 5(0) reproduces semi¬ 
quant it at ively the experimental findings in Fig. 1(c). 

We conclude that the in-plane resistivity anisotropy 
found in the paramagnetic state of iron pnictides depends 
sensitively on the nature of the magnetic phase transition 
and a strong elasto-nematic coupling between the uniax¬ 
ial pressure induced lattice distortion and the electronic 
nematic susceptibility. We caution that while the resis¬ 
tivity anisotropy A p and its dependence on the shear 
strain can be successfully used to extract the quantity 
proportional to the nematic susceptibilit}^, care should 
be taken when equating A p with the nematic order pa¬ 
rameter itself. In particular, the non-monotonic temper¬ 
ature dependence of Ap and its sensitivity to the nature 
of the magnetic phase transition remain relatively little 
explored and deserve further experimental and theoreti¬ 
cal studies. 
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V. APPENDIX 
A. Sample Information 

The iron pnictide single crystals used in present study 
were prepared by self-flux method^. The samples 
have been characterized by resistivity, magnetization, 
and neutron scattering measurements. Figures 5 and 
6 show the basic characterizations of BaFei. 97 Nio.o 3 As 2 


and SrFei. 97 Ni 0 .o 3 As 2 samples, respectively. The basic 
characteristics of the BaFe 2 As 2 and BaFei. 8 sNio.i 2 As 2 
samples can be found elsewherd^ 53 . 


80 90 100 110 120 130 140 



FIG. 5: (Color online) (a) Temperature dependence of 

the magnetic susceptibility under ZFC and FC cases for 
BaFei. 97 Nio.o 3 As 2 - (b) Temperature derivative of the mag¬ 
netic susceptibility, showing clearly the structural and mag¬ 
netic phase transitions, (c) Temperature dependence of the 
Q = (1,0,1) magnetic (green diamonds) and (2,—2,0) nu¬ 
clear (red squares) Bragg peak intensity^*. The extinction 
release of the (2,—2,0) Bragg reflection is sensitive to the 
change of structural distortion and used to determine the T s . 
Combining the results in (a)-(c), the Tn and T s are deter¬ 
mined as shown in the green and red vertical dashed lines, 
respectively. 

Figure 5(a) shows temperature dependence of the zero 
field cooled (ZFC) and field cooled (FC) magnetic sus¬ 
ceptibility x- Figure 5(b) is temperature derivative of 
X, dx/dT. Figure 5(c) shows temperature dependence of 
the magnetic (1,0,1) Bragg peak (green diamonds) and 
(2, —2, 0) nuclear Bragg intensity (red squares). These re¬ 
sults establish T/v (green dashed line) and T s (red dashed 
line) of BaFei. 97 Nio.o 3 As 2 . The T/v was determined as 
109 zb 2 K from magnetic order parameter of the (1,0,1) 
magnetic Bragg peak [Fig. 5(c)] and temperature de¬ 
pendent magnetization measurements [Fig. 5(a)]. The 
structural transition temperature T s is estimated from a 
feature shown in the magnetization, [Fig. 5(b)], and the 
neutron extinction releas e of the (2, —2, 0) nuclear Bragg 
peak intensity [Fig. 5(c) | 41 * 65 ^. 

Figure 6(a) shows temperature dependent resistivity 
R and its derivative dR/dT for SrFei. 97 Ni 0 .o 3 As 2 . The 
dip of the dR/dT at T = 175 K indicates the con¬ 
comitant structural and magnetic transition, different 
from the two features for the separated T/v and T s in 
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FIG. 6: (Color online) (a) Temperature dependent resistiv¬ 
ity and its temperature derivative for SrFei. 97 Nio.o 3 As 2 . (b) 
Temperature dependence of the (4,0,0) nuclear Bragg peak 
intensity measured for warming and cooling. The dramatic 
increase of the (4,0,0) Bra gg r eflection signals the change of 
structural distortion near T s 41 . 


BaFei. 97 Ni 0 .o 3 As 2 . Figure 6(b) is temperature depen¬ 
dence of the (4, 0, 0) nuclear Bragg peak intensity. The 
dramatic increase of the peak intensity also signals the 
structural transition. The observed intensity hysteresis 
is consistent with first order nature of the structural and 
magnetic phase transition. 

Large single crystals were selected and cut into rectan¬ 
gular shapes along the orthorhombic [1,0,0] and [0,1,0] 
directions by a high precession wire saw. The well-cut 
samples were placed inside a uniaxial pressure device 
with b axis being the direction of the applied pressure 4 ^. 
The applied uniaxial pressures for the the samples range 
from P ~ 10 MPa to P ~ 20 MPa, as described before. 

In order to measure temperature and pressure depen¬ 
dence of orthorhombic lattice distortion 5 = (a — b)/(a + 
6), the samples were mounted in the [H,K, 0] scattering 
plane, where orthorhombic (4,0,0) and (0,4,0) Bragg 
reflections can be measured. The effects of uniaxial 
pressure on tetragonal-to-orthorhombic structural transi¬ 
tion and uniaxial-pressure induced lattice orthorhombic- 
ity can be probed directly via measuring temperature and 
pressure dependence of the (4,0, 0) and (0,4, 0) reflec¬ 
tions by neutron Larmor diffraction. For magnetic mea¬ 
surements, the BaFei. 97 Ni 0 .o 3 As 2 sample was mounted 
in the [1,1,2] x [1,—1,0] scattering plane, where both 
the (1,0,1) and (0,1,1) magnetic Bragg peaks can be 
reached 41 . 


B. Neutron resonance spin echo measurements 

Neutron spin echo (NSE) technique has been demon¬ 
strated to be an effective method to measure the slow dy¬ 
namics (quasielastic scattering) with an extremely high 
energy resolution (^ 1 fieV or even to rsj 1 neV)P. By 


combining triple axes spectrometer and neutron reso¬ 
nance spin echo (NRSE) technique, the TRISP spec¬ 
trometer at the Forschungsneutronenquelle Heinz Maier- 
Leibnitz (MLZ) is capable of measuring the lifetime of 
excitations with an energy resolution A E ~ 1 /ieV in the 
range of about 1 — 200 /ieV^. 
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FIG. 7: (Color online) (a) Basic principles for neutron 
spin echo with schematics of the experimental setup. The 
magnetic guide field B directions are clearly marked, (b) 
Schematic diagram for the neutron Larmor diffraction mea¬ 
surements. For NRSE, the neutron precession direction in L\ 
is opposite to that of L2. In neutron Larmor diffraction, the 
neutron precession directions are same in L\ and L2. 

Compared with typical neutron scattering experiments 
where S(Q,uj) (P = huj) is usually measured, neutron 
spin echo measures /(Q, tnse) or P(Q, tnse), where P 
is the polarization of the scattered neutrons, which is the 
time Fourier transform of the 5 (Q,cj) and thus provides 
direct information of S(Q,u j) such as energy line-width 
(lifetime) and intensit^ 48 ^ . 

The basic principle of NSE can be understood in a sim¬ 
plified picture as shown in Figure 7(a). We assume neu¬ 
trons polarized along the y direction with a velocity v\ 
enter the first arm of NSE spectrometer with a constant 
magnetic field B [Fig. 7(a)]. The precession angle in 
the first arm (Li) is then < j>i = = r y\B\Li/vi, where 

7 = 2 /j, N /h = 2.916 &Hz /Gauss is the gyromagnetic ratio 
of neutron, L\ is the length of the first neutron guide arm, 
and t is the time for neutron to travel through the first 
arm. After interactions with the sample, some neutrons 
are scattered into different energy with velocity U 2 . In the 
second arm ( 1 / 2 ), the neutron spin will precess along the 
opposite direction, generating —02 = —As¬ 
suming L\ = L2 = L and V 2 = v\ + 5v, Sv « tq, 
the net phase after passing through both field regions 
will be 0 = Since neutron energy transfer is 
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FIG. 8: (Color online) (a) Measurements of the neutron po¬ 
larization in one period of <j) for r = 16.29 ps and 39.25 ps. 
The solid lines are consine fittings of the data using eq. (8). 
(b) P(r) for various temperatures. The solid lines are fittings 
by eq. (10). The brown dashed line shows a drawing of eq. 
(10) with T = 5 yeV for comparison purpose. The corre¬ 
sponding Lorentzian S(lj) are shown as the solid green and 
brown lines in the inset. 


hjj = \m(v\ — v\) ~ mviSv, the net phase can be writ¬ 
ten as 


0 = 


/ huj L L 
\ mvf 


LU = LdTjSfSE 


where t^se is defined as 


( 5 ) 


t~nse = f — -q~ I = 1-863 x 10 16 B(gauss)L(cm)\ 3 (A) 
\ rnv'l J 

( 6 ) 

Note tnse is not a physical time but a quantity deter¬ 
mined by specific parameters of the spectrometer, with 
the dimension of time. 

The polarization along y direction of the scattered neu¬ 
trons can be analyzed and detected [Fig. 7(a)]. The 
average polarization < a y > for neutrons with energy 
transfer hi j is 


by a factor of 4 compared with the NSE with the same 
B and L. The neutrons only precess in bootstrap while 
keep their spin directions in L\ and L 2 . L 2 can be tuned 
by translating the flipper C 4, by which the intensity with 
respect to the position of C 4, I(x c 4 ), can be measured. 
For a fixed r, the measured intensity can be described as 


I{x c 4) 


h 

2 


1 + P COS 


27r 

A(T C 4 


(x c 4 ^c4,o)] 5 


( 8 ) 


where P is the polarization, Jo is the averaged intensity 
of the scattered beam, Ar C 4 is the period of the intensity 
modulation, and x C 4 ? o is the reference position of C 4. 

The measurements of the P(r) for BaFei .97Nio.03As 2 
under P « 15 MPa are summarized in Figure 8 . Figure 
8 (a) shows the intensity modulations for r = 16.29 ps and 
39.25 ps of Q=(1,0,1) at T = 102 K. The polarizations 
are obtained through fitting the data by eq. ( 8 ). The 
fitted P(r) for different temperatures are plotted in Fig¬ 
ure 8 (b). Assuming the possible broadening in energy of 
the magnetic reflections is caused by some slow dynamics 
(quasielastic scattering), the corresponding S(u) can be 
described by a simple Lorentzian: 


S(u) 


1 r 

7T ((J ~ CJ 0 ) 2 + T 2 ’ 


(9) 


where F (F > 0) (Half Width at Half Maximum) is the 
line-width of the quasielastic scattering (cjo = 0). Fol¬ 
lowed by eq. (9), the P(r) should be fitted by the Fourier 
transform of eq. (9): 

P(t) = P 0 exp( - L), (10) 

All the P(r) in our measurements can be well described 
by this exponential decay, as shown in Figure 8 (b). The 
fitted energy line-widths F are less than 1 peV, meaning 
the signal are resolution limited at the measured tem¬ 
peratures. The comparison between S(u) for resolution 
limited (F < 1 peV) and F = 5 peV is shown in the inset 
of Figure 8 (b), as a reference. The temperature depen¬ 
dence of T is shown in Figure 2(b) of the main text, which 
must have values greater than zero. The large error bars 
for the values of F near T/y are due to experimental un¬ 
certainties of P(r) in Fig. 8 (b). 


C. Larmor diffraction measurements 


< G y >=< cos >= J dujS(Q 1 (j)cos ujtnse- (7) 

Thus < G y > is the cosine Fourier transform of 5 (Q,cj) 
for uj and has been shown equal to the intermediate scat¬ 
tering function J(Q,r). Therefore, the tnse dependent 
polarization P(r), that is, /(Q,r), provide direct infor¬ 
mation about 5 '(Q,cj)PS1. 

In the NRSE, the precession fields and spin flippers are 
replaced by four short bootstrap r-f spin flipper coils [Cl- 
C4 in Fig. 7(b)], which can improve the energy resolution 


We now turn to the neutron Larmor diffraction mea¬ 
surements. Larmor diffraction is a neutron Larmor pre¬ 
cession technique capable of measuring lattice spacing 
expansion and spread with a resolution better than 10 -5 
in term of A d/d. It is sensitive to minor change of lat¬ 
tice spacing d but insensitive to sample mosaicity and 
not much affected by slight misalignment 70 . The total 
precession phase ( (j) to t ) dependent polarizations [P{^tot)\ 
can be used to determine the d spread and the split¬ 
ting between multiple d spacings with small differences, 
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FIG. 9: (Color online) Measurement of the Afoot between 
T — 113 K and T = 118.6 K, by which the relative change of 
d spacing can be tracked. 



<j) tot (rad.) (|) tot (rad.) 



FIG. 11: (Color online) (a)-(d) Precession phase (foot) 

dependent polarizations across the structural transition for 
BaFe 2 As 2 . The solid green curves are fits of the P(foot) by 
one or multi Gaussian d spacing distribution models, (e)-(h) 
are the d spacing distributions in d space. The orthorhombic, 
coexisting two different orthorhombic and tetragonal phases 
can be determined for T — 136 to 139 K, respectively. 


neutron velocity v = hki/m , the total Larmor phase fo ot 
can be written as 


4* tot 


2 ujLLm sin Ob 
ttH 


d. 


(ii) 


FIG. 10: (Color online) Larmor diffraction measurements of 
P(r) for BaFei. 97 Nio.o 3 As 2 . (a) Comparison of P(foot) at 
Q = (4, 0, 0) between T — 4 K and 121 K with P « 10 MPa. 
(b) P(foot) for the (4, 0, 0) and (0,4, 0) reflections measured at 
T =s 102.4 K with P w 10 MPa. The inset shows the Gaussian 
distribution of the d spacing, (c) Temperature dependence of 
the P(foot) for the (4,0,0) peak at P = 0. P(foot ) for the 
two d spacings shows clear modulation which can be fitted 
by eq. (14) (solid curves) due to the twinning caused by the 
tetragonal-to-orthorhombic structural transition. 


such as the peak splitting caused by the tetragonal-to- 
orthorhombic structural transitions in iron pnictided^. 

Figure 7(b) is a schematic of Larmor diffraction. The 
spin flippers are tuned to be parallel with the diffracting 
planes and the neutron precession directions in L\ and 
1/2 are the same. Assuming L\ = = A, the total 

neutron precession phase is fo ot = 2 ujlL/v. From the 
Bragg law |Q| = |G| = 2/^sin##, |G| = 27r/d, and the 


Consequently, the variation of the Larmor phase is pro¬ 
portional to the change of the d spacing (caused by ex¬ 
ternal or thermal effect), that is 

^4* tot = 4>tot—j~- ( 12 ) 

For d change induced thermal expansion, the evolution of 
P(4>tot) at different temperatures and the relative change 
of the foot (A foot) can be obtained by fitting the in¬ 
tensity modulations I (foot, o + A fo using eq. (8), where 
the modulations are measured by scanning x C 4 near x C 4 ; o 
[Fig. 7(b) and Fig. 9]. The Afo ot between different tem¬ 
peratures (or pressures/Qs) can be used to determine the 
evolution of the lattice spacings. To facilitate data analy¬ 
sis, A fo ot between two neighbouring conditions should be 
kept within 27r. Figure 9 is an example of lattice thermal 
expansion at Q = (4,0,0) of BaFei. 97 Nio.o 3 As 2 . The dif¬ 
ference in x C 4 (SL) between T = 113 K and 118.6 K, SL , 
can be converted to lattice expansions according to eq. 
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(11, 12). In present measurements, 1 mm is equivalent 
to ~ 1 x 10 -3 in A d/d, with resolution ~ lx 10 -5 . Note 
lattice expansion measurements is only valid for single 
d spacing at one Q. Systems showing coexisting multi¬ 
ple d spacings around the same Q with small differences 
between them cannot be measured by this method. 

Besides the lattice expansion measurements, the 
P^tot) =< cosA0(0tot) > in a wide range of c/) tot 
has been demonstrated to be the Fourier tran sform of 
the lattice spacing distribution [/(Ad/d) ) 47 | 7Q | . For a 
single-Gaussian distribution of d with FWHM=£^vv 5 the 
P(<t>tot) can be derived as 

P(0 tot) =P o ex (13) 

where the FWHM represents the magnitude of the d 
spread. It is usually expressed in term of A d/d. The data 
shown in the Figure 3 of the main text are temperature 
dependence of the FWHM (lattice spacing spread). All 
of our P((/> to t) fo r single d spacing are well described by 
this model, resulting in a Gaussian distribution of the d 
values. Figure 10 shows the P((j)tot ) of BaFei .97Nio.03As 2 
under P « 10 MPa and their fittings by eqs. (13),(14). A 
clear difference between T = 4 K and T = 121 K in Fig¬ 
ure 10(a) indicates different FWHM of the d spread. Fig¬ 
ure 10 (b) compares P((j)t 0 t) for the (4,0,0) and (0,4,0) 
reflections at T = 102.4 K. The corresponding FWHMs 
of the d distributions are also shown as an inset. Their 
differences suggest that the d spread along the pressure- 
applied orientation is much larger. 

Figure 10 (c) shows temperature dependence of P{(j) to t ) 
for an un-pressured BaFei. 97 Nio.osAs 2 sample. Its evo¬ 
lution at high temperature (T > 109 K) indicates the 
broadening of the d spread. Below 109 K, clear modu¬ 
lations are seen in P(cj) tot ). This is caused by the peak 
splitting of the orthorhombic (4, 0, 0) and (0,4, 0) reflec¬ 
tions in a twinned sample. For systems showing two or 
more coexisting d spacings around some Q, their d spac¬ 
ing distribution functions are superposition of multiple 
Gaussian distributions. In this case, the interference be¬ 
tween different d spacings will appear and can be used to 
identify the specific values and spread of the involved d 
spacings. 

For the peak splitting (two d spacings) in 
BaFei. 97 Ni 0 . 03 As 2 , interference between scattered 
neutrons from d\ and d 2 gives rise to the modulating 
polarization 


P{ 4 >tot) = a 2 + (1 - a ) 2 + 2 a(l - a)cos( 0 tot Ae), 

(14) 

where 

a -%■?>■ < is > 
here we assume both d spacings have the same Spw- a 
and (1 — a) denote the populations of the d\ and d 2 . 


di - d 2 
(d\ + d 2 )/2 ’ 


(16) 
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FIG. 12: (Color online) Temperature dependent orthorhom¬ 
bic lattice distortions for BaFe 2 As 2 and SrFei. 97 Nio.o 3 As 2 . 
The open green diamond in (a) marks the temperature range 
showing four d spacings. The open red diamonds in (b) show 
the persistence of the tetragonal phase into the orthorhombic 
phase, indicative of a first-order structural transition, consis¬ 
tent with previous reports. The vertical blue dashed lines 
mark the structural transitions. 


is the lattice distortion. The definition of A^ is simi¬ 
lar with the orthorhombic lattice distortion S = (a — 
b)/(a + b ) in iron pnictides, with A^ = 2(P. The P(cj) tot ) 
at T — 107.5 K and 105.2 K in Figure 10(c) are well 
described by eq. (14). The fitted lattice distortions and 
d spreads are shown in Figure 3 and 4 of the main text. 
The resolution in determining As here depends on the 
range of (j) to t and the d spread of the sample since the 
dips of the polarization is critical for fitting As. The 
resolution of As for two d spacings is 7 x 10 -4 in the 
present work. Thus the possible distortions at tempera¬ 
tures slightly lower than T s = 114 K [such as the 109.5 K 
data shown in Figure 10 (c)] in BaFei. 97 Nio.o 3 As 2 cannot 
be distinguished from the broadening of the d spread. 


Figure [TT] are P((f>tot) f° r temperatures across the 
structural transition of the BaFe 2 As 2 sample. P( 4 >tot) 
in Figure [Tl|a) is a beating pattern caused by interfer¬ 
ence between two d spacings below T s , similar with that 
shown in Figure 10(c). The corresponding d spacing dis¬ 
tributions are shown in Figure 11 [e). The orthorhombic 


distortions can be determined as 6 = (d 2 — di)/(d 2 + di). 
Upon warming the sample to T = 137 K, a tempera¬ 
ture slightly lower than T s , a more complicated pattern 
[Fig. [Tljb)] indicates the coexistence of four d spac¬ 
ings [Fig. pT[f)]. This is consistent with the coexisting 
orthorhombic antiferromagnetic (di) and orthorhombic 
paramagnetic (S 2 ) phases revealed by high resolution X- 
ray diffraction measurements 5 . The four-d spacing model 
of P(t) can be derived analytically (not shown here) and 
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fit the data very well. Figure [TlJ c) and (g) are results for 
T ~ T s , where the 82 is indistinguishable and only one 
broad d spread can be fitted. Here, the orthorhombic an¬ 
tiferromagnetic phase ((Si) is about to disappear, suggest¬ 
ing this temperature is near T s . For temperature higher 
than T s in Figure [TT|[d) and (h), only one d spacing is 
observed, indicating the system enters into the paramag¬ 
netic tetragonal phase. In Figure [llj a,b), the magnitude 
of the lattice distortions determines the beating periods 
(overall line shape) and the relative populations of differ¬ 
ent d spacings control whether the polarization can reach 
zero at dips. In the present study, the overall line shapes 
of all P(r) are well fitted by specific multiple (2 — 4) d 
spacing models [green curves in Figure [TT|(a)-(d)] , indi¬ 
cating that the lattice distortions are well determined. 

The orthorhombic lattice distortions for BaFe 2 As 2 and 
SrFei. 97 Ni 0 .o 3 As 2 obtained from Larmor diffraction mea¬ 
surements of Pifaot) are shown in Figure 12 These re¬ 
sults are consistent with previous results measured by 
X-ray diffractiorf 5 . The error bars in Figures 3 and 4 of 
the main text are fitting errors of the raw data at different 
temperatures according to formulas discussed above. 


D. d spread anisotropy between a and 6 

Another interesting discovery is the doping dependent 
d spread anisotropy under uniaxial pressure. The sam¬ 
ples shown in Figure 13 exhibit similar temperature de¬ 
pendence of the d spread between (4,0,0) and (0,4,0), 
suggesting that the difference of d spread between a and 
b is trivial. However, we note that the FWHM of (0,4, 0), 
along the uniaxial pressure direction, is much larger than 
a in underdoped samples [Fig. 13(a)-13(c)]. This may be 
attributed to an inhomogeneous distribution of the pres¬ 
sure induced strain field. However, we find very small 
differences in d spread between (4,0,0) and (0,4,0) in 
the overdoped BaFei.ssNio. 12 AS 2 [Fig. 13(d)], suggesting 
the d spread anisotropy between a and b is non-trivial 
and may be associated with antiferromagnetic/structural 
instability or even nematic susceptibility in underdoped 
samples. 


E. The lattice distortions and Young’s modulus 



FIG. 13: (Color online) Temperature and doping dependent 
d spread of BaFe 2 -ccNi x As 2 and SrFei.97Nio.o3As2 measured 
with finite uniaxial pressure. Tn and T s are marked as blue 
(and green) vertical dashed lines. The red dashed line marks 
the temperature with the maximum of the FWHM. 


affect temperature dependence of the pressure-induced 
FHWM of A d/d and its comparison with other iron pnic- 
tides, thus will not alter the conclusions of our experi¬ 
ments. 


The Young’s modulus Y along the 6-axis (^ Cqq) can 
be estimated by Y — P/5, where 5 is pressure induced 
lattice distortion. At ~ 250K, the Y for BaFe 2 As 2 , 
BaFei. 97 Nio.o 3 As 2 and BaFei. 8 sNio.i 2 As 2 estimated from 
our neutron Larmor diffraction experiments are ^ 50 
GPa, ^ 50 GPa and ~ 100 GPa, respectively. Com¬ 
pared with the shear modulus Cqq obtained by ultra¬ 
sound spectroscopy 5 ^, the estimated Y for x = 0 and 
x = 0.03 are ~ 30% larger. These differences are mainly 
caused by the errors in our estimation of the applied pres¬ 
sure P through measuring compressed spring distances 
and estimated spring constant 41 . However, they will not 


F. Landau theory and effect of magnetism on 
nematicity and strain 

In order to understand the distinct behavior of the 
observed lattice distortion in SrFei. 97NI0.03 As 2 compared 
to BaFei. 97 Ni 0 . 03 As 2 [see Fig. 1(c) in the main text], we 
write down the Landau free energy incorporating the 
electronic nematic order parameter </?, coupled magne- 
toelestically to the lattice distortion 5 oc e 6 , as well as to 
the antiferromagnetic order parameter M: 

F = F 0 +Tof[<p,S]+T 0 f[M,'p], 


(17) 
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where 


f\ A1 - ffl T ~ r 0 2 , 

^ 2 T 0 ^ + 4T 0 


4 , ^66,0 ^2 ^ ^ -P£ 

r 2T 0 T 0 r T 0 

(18) 

Here we chose to normalize the free energy by the Curie- 
Weiss temperature To associated with the quadratic p 2 
term (if Tq is negative, it is replaced by |Tq|). Treating 
the electronic nematic order parameter Lp as a dimen¬ 
sionless variable, this has an advantage that all the co¬ 
efficients in the free energy are dimensionless (here we 
choose, without loss of generality, a = B = 1). The re¬ 
maining Landau expansion parameters can be fixed from 
the experiment. Indeed, it is convenient to express the 
external uniaxial stress P in terms of the dimensionless 
stress variable a = P/Cq 6 ,o- Then, the last three terms 
in Eq. (18) can be written as follows: 


A 

W\ 


c 66 , 0 \ ft 2 


— < 5 <t ) — 5<p 


A ) V 2 

Minimizing the free energy with respect to S, we find 

A 


(19) 


5 = 


^ 66,0 


p ± cr. 


( 20 ) 


From the minimization with respect to p, it is easy to 
obtain 


dp A 

d 5 = a(T - T 0 ) ± Wcj ) 2 


( 21 ) 


and now the shear modulus Cqq = dP/dS = d 2 F/dS 2 
becomes 


Cqq Cqq o - A-y = C 66,0 
do 


X 2 


i(T — Tq) ± 3T </> 2 5 


( 22 ) 


in other words the elastic modulus gets renormalized 
from its bare value Cq 6 ,o by virtue of the elasto- 
nematic coupling A. Equivalently, it follows from the 
above equation that the inverse nematic susceptibility 
X ” 1 = d 2 F/dp 2 \p^o also gets renormalized from its bare 
value X p = a(T -T 0 ): 

Xp - X? ~ ^ * a(T - Tf w ), (23) 


where T^ w = Tq ± - is the renormalized Curie- 

Weiss temperature. One can now cast Eq. ( 22 ) above 
the transition temperature T s = T^ w into the form 


d5 _ 1 1 T — Tq 

dP = Cie~ C^T-T?w' 


(T > Tf 1w ) (24) 


which is Eq. (3) in the main text. We now use this 
Eq. ( |24| ) to fit the data for the pressure-induced distortion 
S(P) — (5(0) ~ P(dS/dP) [Fig. 1(c) in the main text], 
extracting the values of Tq and T^ w which we quote in 
Table [I] for BaFe 2 As 2 and SrFei. 97 Nio.o 3 As 2 . 


We now turn to the question of the strength of the 
elasto-nematic coupling constant A. The unknown di¬ 
mensionless parameter r = Cqq$/\ in Eq. (20) can be 
fixed from the ratio p/5 in zero external stress (cr = 
0). Substituting the typical value of 5 ~ 3 x 10 -3 in 
SrFei. 97 Nio.o 3 As 2 and BaFei. 97 Nio.o 3 As 2 and choosing the 
value of the nematic order parameter p = 1 deep in¬ 
side the nematic phase for convenience, we find r ~ 300. 
As for the value of A = A/|T 0 | itself, it can also be 
fixed from the experiment since A enters in Eq. (23) 
to renormalize the Curie-Weiss temperature. Therefore, 


~ / rpC XV rji \ 

one obtains A = a y s ^ 0 ) r. Plugging in the values of 

T^ w and To from our fittings of the lattice distortions 
(Table |I|, we obtain A « 80 for SrFei. 97 Ni 0i 03As 2 and 
A « 490 for BaFe 2 As 2 , in other words the effective 
electron-lattice coupling is about ~ 6 times weaker in 
SrFei. 97 NI 0 . 03 As 2 compared to the BaFe 2 - cc T c As 2 com¬ 
pounds. For BaFei .97Nio. 03 As 2 , the quality of our data 
was insufficient to accurately determine the bare Curie- 
Weiss temperature Tq (we were only able to determine 
T^ w = 88.5 ± 1.0 K). However, from the estimated 
T^ w — To ~ 50 K by the elastic measurement^, we can 
deduce the approximate value of the coupling constant 
A ~ 390, similar in magnitude to undoped BaFe 2 As 2 - 


TABLE I: Curie-Weiss fitting parameters of the pressure- 
induced lattice distortions in the tetragonal state, see Fig. 
1 (c) in the main text. 


Sample 

T° w { K) 

To(K) 

T^ w — To(K) 

BaFe 2 As 2 

134.9 ± 0.3 

51.3 ±8.7 

83.5 ±8.7 

SrFei. 97 Nio. 03 As2 

170.9 ± 1 

135.7 ±6 

35.2 ±6.1 


We now turn our attention to the magneto-nematic 
coupling. On symmetry grounds, nematic order parame¬ 
ter must couple to M 2 (since magnetization breaks time- 
reversal symmetry, and p does not). This can be shown 
explicitly by considering the magnetization M a,b on the 
two sublattices composed of the next-nearest neighbor 
sites of the square lattice, in which case the nematic or¬ 
der parameter couples linearly to the product (M^-Mg) 
|25l28l36l7l| . Note that this conclusion holds indepen¬ 
dently of whether the microscopic origin of nematic- 
ity is purely magnetic 25 28 71 or 01 . blta j3ll33|:M72|73 | The 
magnetic phase transition itself may be intrinsically sec¬ 
ond order, as in BaFe 2 - cc T E As 2 compounds, or first or¬ 
der, as in SrFei. 97 Ni 0 . 03 As 2 , Bai_ a ,(K,Na) a ,Fe 2 As 2 , and 
Cai_ a Xa ;E Fe 2 As 2 . Below we consider both possibilities: 

fi[M,<p] = |(T - T]v,o)M 2 - V -M a + |M 6 - W M 2 (25) 
h[M, <p\ = |(T - Tiv,o)M 2 + V -M a - wM 2 (26) 

Since we are after the qualitative consequences of the 
magneto-nematic coupling, the precise values of the Lan¬ 
dau coefficients are not essential (we take u = v = w = 1 
and fi = 0.1 for concreteness). 
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FIG. 14: (Color online) Change in the lattice distortion as a 
function of temperature, calculated from the Landau theory 
assuming either 1st order (Eq. 25) or 2nd order (Eq. 26) mag¬ 
netic phase transition coupled to the nematic order parameter 
p. The blue dashed curve is shifted to the left by AT = 0.2T o 
for clarity. These results should be compared to the neutron 
data in Fig. lc in the main text. 


FIG. 15: (Color online) Calculated temperature dependence 
of the nematic order parameter p coupled to either first or 
second order magnetic phase transition. Note that in both 
cases, the bare p 4 term is positive, however the nematic tran¬ 
sition is rendered first order (solid curve) when coupled to the 
first order magnetic order parameter in Eq. (25). 


Having introduced the Landau formalism above, we 
now study the effect of the applied external stress P on 
the behavior of the lattice distortion. The calculated 
temperature dependence of S(P) — 5(0) is shown in fig¬ 
ure [14] for the realistic strain P — 20 MPa and is shown 
to depend crucially on the nature of the magnetic phase 
transition. Indeed, the only difference between the two 
curves is the sign in front of the quartic M 4 terms in 
Eqs. (25) and (26), while all the other Landau expan¬ 


sion parameters are kept the same (the two curves are 
offset horizontally for clarity). Note that for small P, 
S(P) — 5(0) « (d5/dP)P is proportional to the nematic 
susceptibility, which is expected to diverge at T^ w ac¬ 
cording to Eq. (|24|). Both curves in Fig. 14 exhibit an 


enhanced nematic susceptibility on approaching T ^ w , as 
expected. The main difference is the shape of the curve 
on approaching the transition, which has a distinct asym¬ 
metric “lambda” shape in the case of the second-order 
magnetic transition, and resembles closely the experi¬ 
mentally measured 5(P) — 5(0) for BaFei .97Nio.o3As 2 in 
Fig. lc (see main text). By contrast, the Neel transition 
is first order in SrFei. 97 Nio.o 3 As 2 , and the experimental 
behavior in Fig. 1(c) is close to the calculated sharp in¬ 
crease seen in our model (solid line in Fig. 14). Therefore, 


the Landau free energy results corroborate our conclusion 
that the nature of the magnetic transition is crucial to 
the observed temperature dependence of the lattice dis¬ 
tortion. 

We note in passing that for sufficiently strong coupling 
constant /i, the magnetic transition becomes weakly first- 
order even if the intrinsic free energy has positive M 4 
term in Eq. (26). This is likely the explanation for the 


observed change of the nature of the magnetic transi¬ 
tion from weakly first order in BaFe 2 As 2 to second order 
upon Co doping^ 7 I3. However for the values of the cou¬ 
pling constants in Fig. [l4j this effect is imperceptible and 


the main difference between the two curves is due to the 
different intrinsic nature of the magnetic phase transition 
depending on the sign of the quartic term in Eqs. (25) 


and ( [26] ) . We have verified that for the significantly larger 
values of the coupling constant (fi > 0.4 in Eq. [26]), it is 
indeed possible to obtain the shape similar to the dashed 
line in Fig. [14] because the magnetic transition becomes 
effectively first order. In either case, our conclusions re¬ 
main intact. 


G. Interpretation of the resistivity anisotropy 

The resistivity anisotropy A p = (p a — Pb)/{pa~\~ Pb) has 
been widely used as a proxy for the electronic nematic or¬ 
der parameter in the iron pnictides 12 20 . However in some 
compounds, in particular in SrFei. 97 Nio.o 3 As 2 , the resis¬ 
tivity anisotropy is vanishingly small immediately above 
T/v [see Fig. 1(b) in the main text], whereas it is much 
larger in BaFe 2 - a; T a; As 2 . This is puzzling because the 
lattice distortion is comparable in both cases [Fig. 1(c)] 
and, according to Eq. (20), one expects the lattice distor¬ 


tion 5 to be proportional to the nematic order parameter. 

To shed more light on this apparent inconsistency, we 
have plotted in Figure [15] the temperature dependence 
of the nematic order parameter p under the uniaxial 
stress P = 20 MPa. The two curves correspond to the 
first- and second-order nature of the magnetic transition, 
respectively, and the Landau parameters were kept the 
same in both cases (except for the sign of the quartic 


term in Eq. 25 and 26). Above the transition temper¬ 


ature, T > T s , the values of p are predictably small, 
but importantly, they are identical in the two cases. In 
fact, the main difference lies in the temperature depen¬ 
dence immediately below T s . From Fig. [15] it would 
appear that in this regime, the nematic order parame- 
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ter should be smaller for the second-order phase tran¬ 
sition, however this is diametrically opposite from the 
comparison between BaFei.g 7 Nio .03 As 2 and SrFe 2 As 2 in 
Fig. 1(b) (see main text), where the magnetic tran¬ 
sition in BaFei .97Nio. 03 As 2 is second order, yet resistiv¬ 
ity anisotropy is much larger. This qualitative observa¬ 
tion can be made sharper by considering Eq. (20), where 
the coupling constant A is estimated from experiment 
to be a factor of ^ 5 larger in BaFei .97Nio. 03As 2 and 
BaFe 2 As 2 compared to SrFei. 97 Nio.o 3 As 2 , whereas the 
elastic modulus Cq 6 ,o ~ 50 GPa is similar in all three 
materials. Then, BaFe 2 As 2 is expected to have at 
least a factor of 5 larger lattice distortion compared to 
SrFei. 97 Nio.o 3 As 2 , assuming that ip is the same in both 
materials. If one now equates the resistivity anisotropy 
with the nematic order para meter </?, as has commonly 
been done in the literature 1 ^QI, then one is forced 
to conclude that must be about 4 times larger in 
BaFe 2 As 2 due to the larger resistivity anisotropy [see Fig. 
1(b)]. Taken together, one would expect the lattice dis¬ 
tortion S to be a factor of ^20 larger in BaFe 2 As 2 and in 
BaFei. 97 Nio. 03 As 2 compared to SrFei. 97 Nio.o 3 As 2 (factor 
of 4 due to larger resistivity anisotropy, times factor of 5 
due to larger A). And yet this clearly contradicts the ex¬ 
perimental evidence in Fig. 1(c), according to which the 
lattice distortion is almost the same in all three materials. 

One possible way out of this dilemma is that the Lan¬ 


dau theory may not be applicable to describe the ne- 
maticity in the pnictides. However, given the excellent 
semi-quantitative agreement that Landau theory pro¬ 
vides for the lattice distortion (Fig. [l4]above) and its well 
documented success describing the elastic shear modulus 
measurements^, such a conclusion is perhaps not well 
justified. Rather, a much more plausible conclusion is 
that resistivity anisotropy is a poor substitute for the 
nematic order parameter. While it is plausible that the 
two quantities are proportional to each other, as follows 
from the nematic susceptibility measurements^!, the co¬ 
efficient of proportionality need not be constant and can 
have a strong temperature dependence (and likely mate¬ 
rial dependence), as suggested recently by Tanatar et al. 
in the recent study on FeSd^. This material displays a 
non-monotonic temperature dependence of the resistivity 
anisotropy with a peak below T s , qualitatively similar to 
BaFe 2 As 2 . Further theoretical and experimental stud¬ 
ies are necessary to elucidate the precise relationship of 
the resistivity anisotropy and the nematic order param¬ 
eter in the iron pnictides and chalcogenides. Direct mi¬ 
croscopic measurements of the nematic order parameter, 
for instance using the angle-resolved photoemission spec¬ 
troscopy (ARPES) to probe the orbital splitting, com¬ 
bined with the uniaxial pressure measurements, would 
be desirable. 
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